Coherence dynamics of kicked Bose-Hubbard dimers: 
Interferometric signatures of chaos 
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We study the coherence dynamics of a kicked two-mode Bose-Hubbard model starting with an 
arbitrary coherent spin preparation. For preparations in the chaotic regions of phase-space we find 
a generic behavior with Flouquet participation numbers that scale as the entire iV-particle Hilbert 
space, leading to a rapid loss of single particle coherence. However, the chaotic behavior is not 
uniform throughout the chaotic sea, and unique statistics is found for preparations at the vicinity 
of hyperbolic points that are embedded in it. This is contrasted with the low log(iV) participation 
that is responsible for the revivals at the vicinity of isolated hyperbolic instabilities. 



One-particle coherence is the most distinct hallmark of 
Bose-Einstein condensation. The observation of matter- 
wave interference fringes served as an unequivocal proof 
for the generation of dilute gas Bose-Einstein condensates 
(BEC) [2] and the loss of their visibility was used as a 
sensitive probe for atom number squeezing j3] and the 
superfluid to Mott-insulator quantum phase transition 
[4] in BECs confined in periodic optical lattices. Coher- 
ent mean-field Josephson dynamics was demonstrated in 
double- well condensates [5] and serves as a starting point 
for the construction of sub-shotnoise atom interferome- 
ters with Gaussian squeezed states [U [7] . 

The single-particle coherence physics of two coupled 
BECs is captured by the two mode Bose-Hubbard Hamil- 
tonian (BHH) model. In this work we study the coher- 
ence dynamics when the hopping term is periodically 
modulated in time so as to produce a chaotic classical 
limit [8]. Specifically we consider a sequence of kicks, 
hence the BHH is formally identical to that of a kicked 
top ISHS]: 




where j x =(a 1 a2+a 2 ai)/2, and jy = (a\a2—a\ai) / (2i), 
and J z — {fii~fi2)/2 are defined in terms of annihilation 
and creation operators for a particle in mode i = 1,2. 
Number conservation hi + fi2 = N, with hi = a\a, im- 
plies angular momentum conservation with j = N/2. We 
note that the same dynamics can be realized by modulat- 
ing the interaction strength U via a Feshbach resonance, 
keeping the hopping term constant [5]. Either way the 
evolution can be obtained by subsequent applications of 
the Floquet operator F = exp(iTKJ x ) exp(—iTUJ 2 ). 

The integrable limit (T — > 0) of Eq.(flJ features a sin- 
gle dimensionless interaction parameter u — NU/K that 
characterizes the dynamics. For u>l a separatrix ap- 
pears in its classical phase-space [TU1 [TT] , and Eq.Q be- 
comes a variation of a pendulum, with the characteristic 
Josephson frequency uj ps (NUK) 1 / 2 [12]. If the period 
T is finite, there is an additional dimensionless chaoticity 
parameter -Kchirikov = T 2 NUK = (ujjT) 2 that corre- 



sponds to the kicking strength as defined in the standard 
map. 

We study the time evolution of the one-particle coher- 
ence S(t) = \(J) t \/j, starting from coherent spin state 
preparations that can be obtained by rotation of the \ j, j) 
state: 

\9,ip) = exp(-itpj z ) exp(-i6J y ) (2) 

In these initial states all particles occupy a single super- 
position of the two modes, with a normalized population 
imbalance S z — cos(0), and a relative phase ip, so that 
S(0) = 1. Experimentally, such states can be prepared 
via a two step process from the coherent ground state 
when u <C 1, in which 9 is set by a coupling pulse and ip 
by a bias pulse [TT] . 

Evolving such preparations, the integrable T — > 
model, like the Jaynes-Cummings model of quantum op- 
tics [13] , exhibits a series of collapses and revivals of its 
single-particle coherence [121 IT4HT6] . These recurrences 
are manifest in the Rabi- Josephson population oscilla- 
tions, as well as in the average fringe visibility, when the 
two condensates are released and allowed to interfere. 

In previous work we have studied the time evolution of 
S(t) in the T — > model, starting with various coherent 
spin preparations [TU [T51 [TB]. The dynamics is deter- 
mined by the expansion \9, ip) = ^2 v \v)(v\9,Lp) of any 
given coherent preparation in terms of the BHH eigen- 
states \v) pertinent to the specified interaction parameter 
u. The effective number of eigenstates that contribute to 
the dynamics is evaluated by the participation number 
M = l/YlvPv> wnere Vv — \k v \9, <p)\ 2 - For example, in 
the strong interaction regime u > N 2 , the eigenstates are 
merely Fock number states, leading to M ~ VN, result- 
ing in the collapse of coherence on a timescale (U vN)~ x , 
and its revival at t = (UN) -1 [17]. Long time dynamics 
in this case reveals the dependence of M on N, corre- 
sponding to effective higher-order interactions |18j . 

The dynamics are far more complex in the Josephson 
interaction regime 1 < u < N 2 , due to the coexistence 
of nearly-linear and highly nonlinear phase-space regions 
[T2T [15] IT6] . Consequently the dependence of M on the 
N is rich and non-universal. For example, the coherent 



2 




jt/2 n 3it/2 2it jt/2 it 3jt/2 2it 



FIG. 1: (color online) Classical phase-space stroboscopic plots 
(top), time-averaged one-particle coherence S (middle), and 
participation number M (bottom) , for all spin coherent prepa- 
rations | (9, ip). The parameters here and later, unless specified 
otherwise, are T = 1, and K — tt/2. The number of particles 
is TV = 200, while u = 5/tt (left), and u = 8/vr (right). 

preparations \tt/2, 0) and \tt/2,tt} are both characterized 
by equal-population, but with very different participa- 
tion numbers, M ss y/u and M sa ^/u \og(N/u) respec- 
tively. Consequently one observes a substantial differ- 
ence in their coherence dynamics [IS]. These differences 
reflect the elliptic versus hyperbolic nature of the classi- 
cal mean-field dynamics at the vicinity of the two points. 
Moreover, for other coherent preparations, with the same 
energy as that of \tt/2, 7r), we find a much larger partic- 
ipation number M ~ V^V [12]. Thus, a purely semi- 
classical picture, going beyond mean-field by accounting 
for the deformation of the initial Gaussian according to 
the classical motion, does not yet provide a satisfactory 
description of the dynamics. 

Outline.— In this work we contrast the previously 
studied intergable case (T — > 0) with the chaotic case (fi- 
nite T), first introduced in the BHH context in Ref. [8]. 
The time-evolution of the single-particle coherence is 
studied, starting from all possible coherent spin prepara- 



tions. We find that coherent states in the chaotic regions 
of phase-space have participation numbers of the order of 
the entire Hilbert space dimension (M ~ N), as opposed 
to the traditional boundary M ~ y/N that holds for the 
integrable system. This leads to rapid loss of one-particle 
coherence, and practically prevents the collapse and re- 
vival near elliptical or hyperbolic points. Furthermore, 
we observe different types of chaotic behavior: Coher- 
ent preparations located on hyperbolic points within the 
chaotic sea have a significantly lower participation num- 
ber compared with those that reside inside the sea for 
which M f=a N/2. The latter agrees with random ma- 
trix theory (RMT), while the former can be regarded as 
arising from so-called scars p~9l421j . 

Dynamics.— We study the quantum dynamics in- 
duced by the Hamiltonian ([!]), starting from a spin co- 
herent state preparation ^ with an arbitrary 9, (p. To 
the extent that an initial spin coherent state evolves only 
to other coherent states (without being deformed), the 
dynamics can be described by the mean-field equations, 
where the spin operators are replaced by c-numbers |10j 
with 0(1/N) accuracy and S(t) = 1 identically. Thus, 
classicality in the /V-body system is synonymous to one- 
particle coherence. In the top row of Fig[l]we show stro- 
boscopic plots of the kicked-top classical dynamics for 
two representative values of the interaction parameter u. 
The left panel depicts a mixed phase-space with regular 
islands embedded in a chaotic sea, whereas in the right 
panel the islands shrink and the motion is chaotic almost 
throughout the entire phase-space. 

While mean-field dynamics assumes perfect one- 
particle coherence, the nature of the classical motion is 
interlinked to its loss. In order to study the dynamics for 
the non-integrable kicked BHH, we iterate the quantum 
state with the Floquet operator F, and calculate S(t) 
over a timescale t u)J , where t is the number of iter- 
ations (T = 1). We then evaluate the long time average 
S = (1/i) J2t The value of S for any initial coherent 
preparation is plotted in the middle row of FigJT] From 
comparison with the classical stroboscopic plots it is clear 
that one-particle coherence is lost completely for prepa- 
rations lying in the chaotic regions of phase-space. By 
contrast the coherence is better maintained in the regu- 
lar island regions with near-unity values around elliptic 
fixed points. 

Overview of the integrable case.— We have previ- 
ously studied the time evolution of S(t) for the undriven 
integrable dimer model (TUl HSJ [TB] , and found several 
different types of behavior. The key to the analysis lies 
in expanding the coherent states ([2| in the eigenenergy 
basis for given parameter values, and semiclassically eval- 
uating the participation number M [12j . Thus we have 
identified the following generic dependence of M on the 
particle-number N, depending on the characteristics of 
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FIG. 2: (color online) Husimi distributions after t = 1000 
iterations, starting with the spin coherent state \ir/2,n/2). 
The interaction is u — 5/V and u = 8/n in the left and right 
panels, with N = 20 (upper) and JV = 2000 (lower). 

the classical motion: 

!yfu elliptic fixed point 

y/ulog(N/u) hyperbolic fixed point (3) 
\/N\og(N/u) separatrix edge 

Consequently the dynamics of the one-particle coher- 
ence is completely different in various regions of phase 
space |15| , and in particular strong revivals are observed 
if the preparation starts at the vicinity of a hyper- 
bolic fixed point. 

Scars.— In order to gain similar understanding of 
the dynamics in the mixed phase-space of the kicked- 
top Eq. , we expand each coherent preparation ^ in 
the Floquet basis, meaning that \v) are now re-defined as 
the eigenstates of the Floquet propagator F. The partic- 
ipation number M is then evaluated using the expansion 
coefficients. The results for all coherent preparations are 
shown in the bottom row of Fig.[TJ for two representa- 
tive parameter values. As expected the participation is 
generally far greater for preparations corresponding to 
classical points within the chaotic sea, compared with 
preparations in regular regions. However, not all the 
points within the chaotic sea have the same participa- 
tion number. As seen in the bottom right panel of Fig[T] 
there are evident 'scars' in the vicinity of the hyperbolic 
fixed points, though they are immersed in a chaotic sea. 

Transition to chaos.— Below we turn our attention 
to the coherent preparation \tt/2, tt/2). As shown in Fig[l] 
it resides on a fixed point of the classical dynamics that 



FIG. 3: (color online) The participation number M and the 
time evolution of the coherence S(t) for the initial preparation 
|7r/2,7r/2). Panel (a): from top to bottom, S(t) for u = 
6.4/tt, 6.64/tt, 6.66/tt, 6.8/tt, with N = 1200. Panel (b): The 
time averaged coherence S as a function of u and N. Panel 
(c): from top to bottom, the participation M as a function 
of N for u = 4.8/tt, 6.6/tt, 6.8/tt, 8.8/tt. Panel (d): The 
participation M as a function of it and TV. 

undergoes a transition from being elliptical to hyperbolic 
as u is increased. Linearizing the kicked top map around 
this fixed point (see supplementary material [3j), we ob- 
tain the Lyapunov instability exponent, 

A = i- log (b + - l) , (4) 

b = i^-^J [1— cos(ttu)] + ^— J sin(TTu) + cos(7ru) 

where To — 4 is the fixed point period, i.e. the number 
of kicks required to cycle it to its initial position. From 
Eq.Q it is clear that the fixed point becomes hyperbolic 
at u s» 2.1. Note that for larger values of u there are 
narrow windows around u = 4, 6, 8, corresponding to 
2im kicks, within which elliptic behavior is restored. 

In Figj2] we display the time evolved Husimi distribu- 
tion (9, ip\p(t)\9, ip) with p being the instantaneous den- 
sity matrix after many (t = 1000) kicks, for represen- 
tative values of u and N. From the lower panels it is 
evident that for sufficiently large N there is a transition 
from conservation of coherence over long times when u is 
such that the fixed point is elliptic (panel c), to a rapid 
loss of coherence due to the smearing of the distribu- 
tion throughout the chaotic sea when u is larger and the 
preparation is inside it (panel d). Since 1/N serves as 
an effective Planck constant for the many-body system, 
this classical-quantum correspondence is blurred for low 
particle numbers (top panels), when the "size" of the 
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FIG. 4: (color online) (a) The dependence of M on TV for 
three preparations in the case of the integrable model: at el- 
liptic point |7r/2, 0} (bottom A), at hyperbolic point |-7r/2,7r) 
(middle o), and on the seperatrix edge |0,0.97r) (top V). In 
(b) we contrast the results that are obtained in the case of 
the kicked system, for preparations that reside on the hyper- 
bolic point \n/2,n/2) (blue o) and on a nearby chaotic point 
|vr/3, tt/3> (black □). For both (a,b) u = 9/tt. Panel (c) 
displays the dependence of M on u in the two cases of (b) 
with N = 1000. The dashed and dotted red lines depict the 
RMT estimation M m N/2 and the semiclassical prediction 
of Eq.([5|, repectively. 



preparation is large compared to the Planck cell. 

The resulting one-particle coherence dynamics is 
shown in Fig.[3^i. In agreement with Fig^2] coherence is 
maintained for values of u below the elliptic-to-hyperbolic 
transition, and it is abruptly lost above it. The time- 
averaged coherence S is plotted throughout the (u, N) 
parameter space in Fig[3]3, highlighting the sharp change 
at the classical transition, and the loss of quantum- 
classical correspondence at low N. 

A similar sharp transition is observed for the partici- 
pation number M of this preparation, as seen in Fig[3j;,d. 
For small u the fixed point is elliptic, and hence M has 
no dependence on N, as implied by Eq.([3]). However with 
the onset of chaos, when u is increased, we obtain a differ- 
ent linear dependence of M on N. Recall that N + 1 with 
N = 2j is the dimension of the many-body Hilbert space, 
hence we have here an agreement with the semiclassical 
quantum-ergodic picture: a minimal-Gaussian represent- 
ing the initial coherent state, has a non-zero overlap with 
all the eigenstates that reside in the chaotic sea. 

Wavefunction statistics.— Fig. [3] summarizes the 
entire range of participation number behavior studied so 
far, for the various initial coherent preparations. Figj4^, 
includes the previously studied cases for the integrable 
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FIG. 5: (color online) Wavefunction statistics for the two 
representative chaotic sea preparations: The local density of 
states is plotted for \n/2, tt/2) and for \n/3, tt/3) in panels (a) 
and (b) respectively. Note the different vertical scale. The 
intensity statistics is compared in (c): blue (o) and black 
(□), respectively. The dashed and solid lines are based on 
the RMT and on the Semiclassical theories (see text). The 
parameters are N = 1000 and u — 9/-7T. 



BHH model, whereas Figj4j} shows the linear scaling that 
we observe in the chaotic case. The large M that char- 
acterizes chaotic sea preparations is responsible for the 
observed rapid loss of the one-particle coherence in Fig[2] 
The observed scars in Fig[l] are characterized by a sig- 
nificantly smaller M/N ratio. In Fig.|4j3,c we compare 
the participation ratio at the hyperbolic point (produc- 
ing the scar) with the participation ratio of a nearby fully 
chaotic point. In the chaotic sea the RMT expectation 
M « N/2 (dashed red line in Fig.^) is realized [52]. 
By comparison, in the vicinity of the hyperbolic point 
we expect this result to be semi-classically suppressed as 

nasi], 



M 



E 



i 

cosh(As) 



N 

T 



(5) 



Substituting the Lyapunov instability exponent A from 
Eq. Q with no adjustable parameters, we obtain very 
good agreement between this estimate and the hyperbolic 
point participation ratio (see FigJIp). 

The local density of states W{E V ) = \(v\9,tp)\ 2 pro- 
vides a more detailed statistical information on the par- 
ticipation of the eignstates in a given coherent prepara- 
tion. Here exp[— i2%E v \ are the eigenvalues of the Flo- 
quet operator F. Fig.[sJ compares W{E V ) for the two 
coherent preparations: hyperbolic versus generic chaotic 
point. One observes that the state \tt/2,tt/2) projects 
preferably onto a subset of eigenstates, reflecting that 
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the latter have a significantly larger weight at the fixed 
point ("scarring"). 

To better quantify the eigenvector statistics of the two 
representative preparations, we plot the inverse cumula- 
tive histogram: F(w) is the fraction of eigenstates with 
intensity W(E V ) > w. As expected the RMT statistics 
follows the Porter Thomas F(x) — exp(— x) law, where 
x = Nw. In contrast the statistics in the vicinity of the 
hyperbolic point follow a significantly different functional 
dependence [3U1 [ST], namely F(x) oc a;" 1 / 2 exp(— 72;), 
where 7 is a fitting parameter that is proportional to 
the instability exponent A. 

Summary.— Considering the coherence dynamics of a 
non-integrable kicked-top BHH, one observes that an ini- 
tial spin coherent preparation that reside in the chaotic 
regions of the mixed phase-space, contains 0(1) fraction 
of the quantum eigenstates. This leads to an abrupt ir- 
recoverable loss of single particle coherence, as opposed 
to the collapse and revival dynamics that is obtained for 
the integrable model. Within the chaotic sea we find two 
distinct types of linear dependence, reflecting different 
wavefunction statistics. Namely, the lower participation 
number in the case of a wavepacket that is launched at a 
hyperbolic point reflects the wide distribution of overlaps 
due to scarring. 
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Supplementary material 



Consider the Kicked top map with T = 1, and K = tt/2, and U = w/2 such that u = (2/ir)w. Using the notation 
S = (X, Y, Z) , the classical map is 

X' = Xcos(wZ)-Ysin(wZ) (6) 
Y' = Z (7) 
Z' = -Xs\n{wZ) -Ycos{wZ) (8) 

The fixed point \ir/2, n) is cycled from S = (0, 1, 0), to (0, 0,-1), to (0, —1, 0), to (0, 0, 1), repeatedly, and hence has 
period T = 4. The linearized transformation of (X, Z) involves the matrix 

/ cos(2w) + w cos(ui) sin(ui) — 2w cos(w) 2 — w 2 cos(w) sin(w) + w sm(w) 2 + sin(2ui)\ , . 

I — sin(w)(2 cos(w) + w sm(w) cos(w) 2 + 3w cos(w) sin(w) + (w 2 — 1) sin(w) 2 J ^ ' 

The Lyapunov instability exponents ±A equal the log of eigenvalues of this matrix divided by To- 

Note that the fixed-point undergoes a transition from being elliptical to being hyperbolic at u w 2.1. For larger 
values of u there are narrow windows around u = 4, 6, 8, corresponding to 2im kicks, within which the fixed-point 
is elliptic. This is reflected in the calculation of the semiclassical suppression factor 

oo ^ 



s— — oo 



Together with the RMT statistical factor Frmt = 2, that is implied by the Gaussian Unitary Ensemble, it gives the 
overall suppression: 



M « 1 N (11) 

-Frmt Fsc 



As confirmed in Fig|4j:. 



